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Abstract. - We study the influence of the non- homogeneity of a magnetization fleld on the 
behaviour of interacting electrons in a quantum dot. In particular we investigate the magneto- 
transport properties when the dot is weakly coupled to two ferromagnetic leads. We take into 
account the interactions in the quantum dot non-perturbatively. For a magnetization which varies 
slowly on the scale of the Fermi wave length, the non-homogeneity effect is described by a gauge 
potential that can be treated perturbatively. 



> 

o 
o 

o 

00 

o 



X 



. Introduction. — The role of magnetization and inter- 
action in determining the electronic and transport prop- 
erties of quantum dots has already received considerable 
attention. In particular, many efforts have been devoted to 
the regime vsrhere the dot can be considered as a tvsro-level 
System [1-5] or as a double island [6]. On the other hand, 
r.ecent experimental activities indicate a delicate role of 
the magnetization on the magnetotransport properties of 
a metallic dot. E.g., up to several hundred per cent ballis- 
l;ic magnetoresistance was measured for Ni [7-9], Co [10] 
and Fe nano-islands connected to ferromagnetic leads (cf. 
Refs [11,12] for an overview and further references). Such 
metallic, magnetic nanodots can no longer be reduced to 
a two-level system [13, 14]. 

In addition the magnetization may well not be homoge- 
nous, in particular when the exchange length is on the 
scale of the dot size. Hence, we inspect the scenario of 
a nonhomogcncously magnetized, interacting dot with a 
mean level spacing, (5, which is much smaller than all other 
relevant energy scales. At suitably low temperatures such 
quantum dot systems display the Coulomb blockade ef- 
fect [15]. Specifically, we consider two non-coUinear ferro- 
magnetic leads coupled to the quantum dot and investigate 
the transport properties. Such a set-up is also suited to in- 
vestigate scanning tunneling microscopy measurements on 
islands which have non-homogenous magnetizations [16]. 



As far as we are aware such a system had yet to be studied. 

The quantum dot may be viewed as the domain wall 
region between the ferromagnetic leads. For low-density 
diluted magnetic semiconductor wires we studied a similar 
situation using a Luttinger-liquid and a renormalization 
group approach [17] under the assumption that the car- 
riers Fermi-wavelength is larger than the dot region; the 
effect of which can then be modeled as a point-like, spin- 
dependent scatterer. These conditions are not applicable 
to the metallic case, because of the much smaller Fermi 
wave length the influence of the details of the magnetiza- 
tion profile need to be considered in this case. 

The model. — We consider a quantum dot with a 
spatially non-uniform magnetization M(r), which couples 
to the electron spin density with a strength J. The dot is 
weakly coupled to the ferromagnetic leads that have uni- 
form magnetizations aligned in opposite directions. Due 
to the weak coupling to the leads we consider, to lowest 
order in the coupling, the quantum dot as being closed. 
The Hamiltonian for the dot region in terms of the free 
and interacting contributions reads 

Hq = I Avai{v)i{v)aa{v) -J draicra0.M{r)ap (1) 
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We used units in which ft = 1 and fcs = 1 throughout. 
Repeated indices are summed over. is the creation 

operator for an electron of spin a and with a quantum 
number k which labels the states in the dot. ^ — e — ii 
is the dispersion for the quantum dot measured from the 
chemical potential. This chemical potential includes a gate 
voltage applied to the quantum dot. Eq. (2) describes 
a simplified Coulomb interaction [18, 19] in a quantum 
dot with a charging energy Ec ~ e'^ /2C, C is the self- 
capacitance of the dot. A vector gauge transformation is 
performed to simplify the magnetization term, allowing 
us to use perturbation theory on the resulting introduced 
potential. After this transformation we will have a Zee- 
man splitting term and a spin-dependent spatially vary- 
ing potential, C/a/3(r) [20-23]. Perturbation theory is valid 
provided the magnetization varies slowly compared with 
the Fermi wavelength of the electrons, which is usually the 
case for a metallic dot. The gauge transformation [20-23] 
is 



in the functional integral representation for the Green's 
function with the above Hamiltonian [25], in the Keldysh 
representation [26,27], we have 
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with the action given by 
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where the unitary matrix T(r) is defined such that 

T\r)a{r).n{r)T{r) = cr^ where (4) 
M(r) = n{r)M. (5) 

n(r) is a unit vector and M is the, spatially invariant, size 
of the magnetization. After the transformation our new 
Hamiltonian can be written as 



We define ^ka = Ska-fJ- = ek~ fJ-- JMcr^a and k= {k,a}. 
The contour c is the un-rotated Keldysh contour [14] and 
all times are defined along this contour. 

We can rewrite the effect of the interaction using a 
Hubbard-Stratonovich transformation [14], introducing a 
bosonic field 6: 



P is the inverse temperature. Now 
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with the potential given by 



1 



iSe = iJ2 [ dtMt)[^dt-^k + 0/P]Mt) 

1, J c 



U(r) = _--^[2A(r).9r + ^r-A + A^ir)]. 
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The vector potential is defined, in terms of the transfor- 
mation, as A(r) = Tt(r)apT(r). 

For the case in which the magnetization is translation- 
ally invariant in the x and y plane we set n(r) n[z) and 
parameterize in terms of an angular fimction (p(z): 



Z{e) = / V>i:V)i:e^o = W[l 



(14) 



The average in eq. (13) is defined as 
/ d9e^^ . . . 
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Micromagnetic simulations for a magnetic stripe of a 
length L and a width and a thickness smaller than the 
exchange length deliver the angular profile ^p{z) = tt — 
cos~^ [tanh(z/L)] (cf. Ref. [24] and references therein for 
further details and experiments). Thus eq. (7) takes on 
the form 



This is our starting point for perturbation theory. After 
expanding in U we can perform the averaging over 0, this 
is achieved by transforming the grand canonical ensemble 
into a sum over canonical ensembles. 
Let us first define the contour 
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Firstly we must manipulate the interacting term so that 
our Hamiltonian is quadratic instead of quartic. Working 



In standard diagrammatic perturbation theory [28] we 
have, to first order: 

see eq. (17). 
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Where we have additionally defined 
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We need to calculate terms such as 

(Zkk'(^)e^^^^-^^=-^'').. (20) 
We transform to the sum over canonical ensembles thusly: 
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(We have introduced = E,,N^/2 - nN.) The function 
•^Af(£k,£k') can be rewritten as 
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[l-Fjv(£k)][l-i^w(£k')]- (23) 



^Af(ek) is the canonical iV-particlc distribution function 
for being in any TV-particle state containing the level Ek- 
For large N we find FAr(ek) ~ /(^k), the Fermi distribu- 
tion. Combining all of the above results thus far we can 
write the Green's function to first order: 

see eq. (24). 

iG'f.^{t,t') is the Coulomb blockade result in the presence 
of a Zeeman splitting term [14] and 



(25) 
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The second order terms are calculated in the same manner. 

The density of states. — The standard formula for 
the density of states in terms of the advanced and retarded 
Green's functions is given by 

(^) = 1- trfe [zGf^,,^ (a;) - zG^^,,^ (u)]. (26) 



We find no first order correction in U to the standard 
Coulomb blockade result [15], only second order terms. 
For small distances compared with the size of the dot L, 
|z| <^ L, the variation of the magnetization inside the dot 
is approximately ip{z) ss 7r/2 -I- z/L. This approximation 
still remains very accurate up to |z| ^ L and allows us to 
perform the necessary sums and integrals. Hence we use 
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m is the electron mass. The magnitude of the magnetiza- 
tion is taken to be such that JM > Ec- 

By varying an applied gate voltage to the dot we can 
shift the system through the Coulomb blockade valleys 
and peaks. We write the effect of the gate voltage as 

= Ec{No + i) + (5^. A^o is some (large) integer and 
measures the distance from the degeneracy point, i.e. the 
peak in the Coulomb blockade regime. Sfi = Ec/2 is situ- 
ated at the centre of the Coulomb valley. This pattern re- 
peats periodically in Ec- The density of states is depicted 
in figure 1. Figures 2 and 3 show the corrected density 
of states in comparison with the zcroth order density of 
states for different values of the applied gate voltage. At 
this order in U there is no difference between the density 
of states for spin up and for spin down. 

We note that there is no longer any complete suppres- 
sion of the density of states at any point. The scattering 
opens up additional states on the scale of JM around the 
Fermi level. Below the Fermi energy we note the appear- 
ance of subsidiary structures in the density of states. On 
closer inspection it becomes clear that their origin lies in 
the scattering from states at an energy 2JA/ below the 
Fermi level. 

Spin-dependent current. — For studying the trans- 
port properties we add two quasi-one dimensional mag- 
netic leads attached to the quantum dot by point con- 
tacts. The quantum dot is described by the Hamiltonian 
H — Hq + Hj + Hit as before. We introduce a coupling 
term, Ht^ between the dot and the leads. The current 
through the dot is given by [29,30] 

I ^ ei[H, N] = ei[Hu N] (29) 

El — Hdot + Ht + Hieads (30) 

Ht = ^[tankrfi„ak + Cnkflk'^an]- (31) 
ank 
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Fig. 1: A contour plot of the density of states v{e,5yL) /uq as a 
function of the energy (e/Ec) and the gate voltage 5^./Ec. 



Where a labels the leads (left and right), is the creation 
operator for electrons in the dot and d]^ is the creation 
operator for electrons in lead a; the Vs describe tunneling 
between the dot and the leads. The current reads 



x[iG^k'H-*Gi^k'H]r£'kH, (32) 

^ 27rz^"„ {taa"M't*aa" ,fc V ) > (33) 



where we have assumed t to be independent of the lead 
states. iG^^^,^,{uj) is the dot Green's function coupled 
to the two leads. To lowest order in i^^T" it is the un- 
connected dot Green's function, i/" is the density of 
states for lead a. The distribution function for lead a 
is h"{u}) = 1 — 2/"(a;) whereas h'^{(jj) is the dot distri- 
bution function. Demanding the flow of electrons in and 
out of the quantum dot to be balanced determines the 
steady-state lead distribution function and the current is 
then 
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Fig. 2: The density of states ^{e, Sfi) /vq (dotted) compared 
with the zeroth order form (full), as a function of energy e/Ec- 
The energy, e, is measured from the Fermi energy. The gate 
voltage is tuned here to the Coulomb blockade peak. The in- 
set shows a case intermediate between the Coulomb blockade 
valley and peak. 



We linearize the dispersion relation near the Fermi en- 
ergy and take F to be a constant in /c-space, equivalent 
to assuming all momentum states in the dot are equally 
correlated. Furthermore we assume that there is no spin 
scattering on tunneling and we have two ferromagnetic 
leads of opposite spin orientation. For simplicity we as- 
sume that F^ = F^. 

The linear differential conductance, for small biases, is 
given by 



dl 
dV 



(35) 



V=Vl-Vr=0 
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The results are shown in figures 4 and 5. We plot the 
scaled differential conductance G/Gq where Go = e^Fi/g 
is the conductance for temperatures much larger than the 
energy when charging effects inside the dot will play 
no role. For a metallic system of size L « 10~^m, E,. « 
20meV. The temperature is taken to be « 23K, a tenth 
the size of the charging energy E,.- Also Ec ~ U . 

As can be seen from figure 5 the structure of the second 
order correction to the linear differential conductance, see 
figure 4, is too small to be observed. The main feature 
we note is that the Coulomb blockade valley is weakened, 
there is no longer a complete suppression of current at 
these points. The reason for this is clear if we look at 
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Fig. 3: The density of states i/{e, Sfj,) /vo (dotted) compared 
with the zeroth order form (full) as a function of energy e/Ec- 
The energy, e, is measured from the Fermi energy. The gate 
voltage is tuned in this case to the centre of the Coulomb block- 
ade valley. 
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Fig. 4: The second order correction to the linear differential 
conductance as a function of 6fi/Ec. We plot G^'^^/Go- The 
inset show the magnetotransport against We define 

qM ^ (GTi_G")/G", G^'^ being the differential conductance 
between leads with a spin direction of 5 and 7. 



the density of states in the valley, figure 3. The scatter- 
ing from the potential U introduces states into the pre- 
viously empty region and hence it is always possible for 
electrons to tunnel through the dot. We emphasize that 
the difference between the curves in figure 5 is absolute. 
Thus it is possible to measure an absolute difference in the 
conductance due to the non-homogenous magnetization. 
Additionally, although the structure of the second order 
correction is small, it can still be observed in the magneto- 
current. The inset of figure 4 shows the magneto-current: 
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G^-'- is the current between non-coUinear wires and G^'^ 
is the current between coUinear wires, i.e. the standard 
Coulomb blockade result with a homogenous magnetiza- 
tion in the dot. 

We will also consider the non-linear differential conduc- 
tance. Let us set Vr = and look at 



G(<5Ai, V) 



dl 
dV 



(37) 



Vr=0 



This offers an approximate way of viewing the actual 
structure of the density of states [31] , which becomes exact 
at zero temperature. By tuning the gate voltage and the 
bias voltage we can map out the whole density of states, 
compare figures 1 and 6. This pattern repeats periodically 
with the gate voltage, Sfi. 

Conclusions. — We have investigated the interplay of 
the Coulomb blockade and a non-homogenous magnetiza- 
tion for a quantum dot coupled to non-coUinear ferromag- 
netic wires. We have shown that the scattering from the 
magnetization reduces the effect of the Coulomb blockade. 
By calculating the full profile for the density of states we 



clearly see the origin of this effect. States are created, by 
electrons scattering from the magnetization profile, in the 
previously depleted portion of the density of states. This 
occurs within an energy gap of the order of JM from the 
Fermi level. Additionally we were able to calculate the full 
non-linear differential conductance, which would allow, in 
principle, the full mapping out of the density of states. 
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Fig. 6: The non-linear differential conductance, G/Go, as a 
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